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We study the ground state properties of an anisotropic, quasi-2D version of the quantum (S=l/2) pyrochlore 
antiferromagnet. In the presence of Dzyaloshinsky-Moriya interactions, in addition to the Heisenberg ex- 
changes, it is shown that two types of ordered magnetic states are generally possible: non coplanar "chiral," 
and coplanar antiferromagnetic order. The magnetic moments in all cases are determined by the Dzyaloshinsky- 
Moriya interactions and in this sense the antiferromagnetic order is "weak." 



I. INTRODUCTION AND DESCRIPTION OF THE MODEL 

The pyrochlore antiferromagnet is strongly frustrated and 
the structure of its ground state represents a challenging prob- 
lem in the field of magnetism. It has been argued that, for 
purely Heisenberg interactions, both the classical and quan- 
tum versions of the model are magnetically disordered. '•^•■^ 
It also appears that in the extreme quantum limit (S=l/2), 
lattice-symmetry breaking and spontaneous dimerization take 
place, although the ground state still exhibits macroscopic 
degeneracy.-'''*'^ Thus it is natural to expect that interac- 
tions beyond the Heisenberg exchange, as ■well as inclusion 
of various magnetic anisotropics, can impose their own or- 
der, such as orbital, dimer, or magnetic, or combinations of 
these. Effects due to orbital degeneracy,^ long-range dipolar 
interactions,^ and spin-lattice interactions^ have been studied. 
Ising anisotropics can lead to the formation of "spin ice",^ 
which bears resemblance to the problem of proton disorder in 
ice.'° Various planar"-'- and other anisotropic '-^ versions are 
also of strong theoretical interest, although they do not neces- 
sarily reflect the physics of the full 3D pyrochlore structure. 
Finally, under certain conditions, even more exotic ground 
states have been proposed, such as U( 1 ) spin liquids, *'* or non- 
magnetic chiral states.'^ 

In this work we will discuss certain aspects of the mech- 
anism for magnetic order formation due to the presence of 
Dzyaloshinsky-Moriya (DM) interactions."' For the case of 
classical spins, this mechanism was studied recently'^ (see 
also'^). The quantum case (S=l/2) was discussed'^ within a 
technical scheme, similar to the one used in Ref. 3 for the pure 
Heisenberg case. A puzzling difference between the classical 
and the extreme quantum case is that in the classical case non 
planar, chiral-like (see below) as well as coplanar spin order 
is possible,'^ while in the quantum case only the chiral or- 
der was predicted. ''^ Naturally, the way spin order emerges is 
also quite different in the classical and quantum cases, and 
for S=l/2 the induced antiferromagnetic order is "weak," in 
a sense that the ordered moment is proportional to the DM 
interaction itself. We point out that the role of DM interac- 
tions has been extensively studied only in the context of non- 
frustrated lattices, such as the square lattice, where it typically 
leads to weak ferromagnetic moments (present, for example, 
in the copper oxide compounds.) In frustrated lattices, how- 
ever, the DM interaction effects can be much more profound, 
and are expected to lead to complex types of order ^" 

We consider the anisotropic version of the 3D pyrochlore 
lattice (Fig. 1(a)), where it can be viewed as weakly-coupled 




FIG. 1: (Color online) (a) Pyrochlore lattice (not all tetrahedra 
shown, for clarity.) (b) 2D Pyrochlore "layer," showing the bonds 
in a (001) slice of the pyrochlore lattice, (c) DM interactions on a 
tetrahedron. The DM vectors are represented by blue arrows. 



"layers," defined in Fig. 1(b). Our goal is to investigate the 
possible types of magnetic order arising in this situation. It 
is technically very advantageous to consider the anisotropic 
limit, namely Jj^ <C J' <C J, where Jj^ is the Heisenberg 
exchange between the layers (thin black lines in Fig. 1(a)), 
and J' is the inter-tetrahedral exchange (solid blue lines in 
Fig. l(a,b)). J is the exchange on the "strong" (shaded) tetra- 
hedra, shown as plaquettes in Fig. 1(b). This strong-coupling 
approach is similar to the one used earlier,^'*'''' except that 
now we consider an anisotropic version of it. The quasi 2D 
version allows us also to monitor the low-energy excitation 
spectrum, and determine the conditions under which different 
types of DM-induced order can arise. In order to implement 
this program we need to know the exact excitation spectrum 
on a single tetrahedron, which we calculate below, and then 
consider the lattice version of coupled tetrahedra. The lattice 
shown on Fig. 1(b) is a frustrated one, and we will present 
arguments why the strong-coupling expansion, governed by 
the parameter J' < ,/, should work well in determining the 
ground state structure. In order to avoid cumbersome for- 
mulas we present below results for the strict limit Jj^ = 0, 
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while we have checked that the types of order we find are 
pretty generic (as long as the system is away from the strictly 
isotropic 3D case, which requires different considerations.) 
Notice also that the lattice of Fig. 1(b) is not the same as 
the 2D projection of the 3D pyrochlore, i.e. the checkerboard 
lattice," which has valence-bond soUd order (expected, quite 
generally, to compete with DM-induced magnetic order.) 
The Hamiltonian reads 

H = X]^i.jSi.Sj+X]Di.j.(SixSj), (1) 

ij ij 

where the couplings Jj.j are assumed to be antiferromagnetic 
( Jij > 0) and are distributed as already discussed. All spins 
are S=l/2. The DM vectors Dj j on a tetrahedron are shown 
in Fig. l(c),'^ and have equal magnitude [Dj j| ~ \D\. Two 
patterns are possible: the one shown on the figure, and a pat- 
tern with all DM vectors reversed, Di j ^ Di.j- There is no 
reason to expect that the induced order in the two cases will 
be the same; in fact two different types of order were found in 
the classical version. 

In the rest of the paper, we first calculate the exact spectrum 
on a single tetrahedron (Section II), which is then used to an- 
alyze the spectrum on the lattice and determine the possible 
types of order in Section III. Section IV contains our conclu- 
sions. 




FIG. 2: (Color online) (a) The gaps A = £5 -Bo (solid line), Ah = 
Eb{H) — Eo (dashed line), and the "hopping" parameter t (Eq. (18)), 
solid blue line, as a function of the DM interaction strength, (b) 
The parameter F (Eq. (15)) which determines the magnitude of the 
ordered moments. D is in units of J = 1. 



II. QUANTUM SPINS ON A TETRAHEDRON WITH DM 
INTERACTIONS 



Without DM interactions, the ground state on a tetrahedron 
of spins S=l/2 is a twofold degenerate singlet. We denote the 
two states by |si), |s2). Their explicit definition is given in 
Appendix A. The spectrum above the ground state consists 
of 3 degenerate triplets {\pa),\qa),\ta), a — x,y,z), see 
Appendix A, as well as S=2 states which are irrelevant for our 
purposes. However, the DM interactions break spin-rotational 
symmetry, and thus lead to mixing of singlet and triplet states. 
The ground state is still degenerate, and the two new ground 
states are: 

I*) = a\s2) 

+ ^ {\P.) + \Py) + Ife) - +2a/2|<,)} , 
1$) = a\si) + i/3 - \py) + \q,) + \qy)) . (2) 

From now on we measure all energies in units of J, i.e. we 
set J = 1, and consequently J'/ J J', D/J D, etc. 
The ground state energy, corresponding to the states (2) is 



Eo - -1 
For small D ^ 1, 
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The coefficients a, in the wave-functions (2) are given ex- 
plicitly by the formulas 



" V^Ti' A' ~ 3/2 + Eo- 
It is also useful to have the expansion for small D, 

a^sign{D), (3 ^ ^\D\ 
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(5) 



(6) 



In our convention, the positive sign of D, D > in all the 
above equations corresponds to the pattern shown in Fig. 1(c) 
(called "indirect" in Ref 17), while D < ("direct" case 
from Ref. 17) is the situation when all arrows in Fig. 1(c) are 
reversed. \D\ is the magnitude of the Dij vectors. 

Next, the first excited state (which we call |5)) is 



l^> = ^ {b.) + \Py) + Ife) - \qy) - V2\U)} 



V6 

with the exact energy 
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Eb = -- - V2D. 



(7) 



(8) 



The next excited state is the "triplet" |P) with components: 
1 



IP. 



{-\P=c) + \Py) + \q.) + \qy)} ■ 



\P.,y) = ^{\p,)±\q,)+V2\t,,y)} 



(9) 
(10) 



The corresponding energy is 



1 \/2 
Ep = —D 



(11) 



For _D = all excited states are degenerate and separated 
by an energy gap J = 1 from the ground state. However 
for finite D, the 6-state is the lowest (for D > 0), and the 
variation of its gap A = Et ~ Eq as a function of D is shown 
in Fig. 2(a). In fact the gap vanishes for the (unphysically) 
large value of Z? = a/2; on the other hand the gap increases 
for D < Q. We proceed to explore if this difference can affect 
the type of magnetic order on the lattice (i.e. we study the 
conditions under which b can condense.) 



III. LOW-ENERGY EFFECTIVE LATTICE THEORY: 
CHIRAL VERSUS COPLANAR SPIN ORDER 

A. Low-energy dynamics of interacting spins 

The ground state structure from (2) implies that the spin op- 
erators have non-zero matrix elements of the type (\I'|S,j|$). 
This is due to the presence of triplet states within the ground 
state subspace. Using the explicit form for the various triplet- 
triplet transition matrix elements, summarized in (Al), we ob- 
tain 

S, = r A, 7;^ (12) 

where the vectors are defined as 

Al = (-1,-1,1), A2 = (-1,1,-1), (13) 
A3 = (1,1,1), A4 = (1,-1,-1). 

The site label i = 1,2,3,4 in S;, A;, refers to the number- 
ing of spins on a tetrahedron as displayed in Fig. 1(c). The 
operator is defined as (its bold index i refers to the whole 
tetrahedron) 



(14) 



where $^ are operators that create the two ground states. 
Obviously is the y ("magnetic") component of the pseu- 
dospin operator T = 1/2, defined in such a way that = 
±1/2 label the two ground states (2), i.e. = i(<I>'''$ - 
vl'1^5'). The coefficient F in (12) depends on D through the 
wave-function coefficients. 



af3 - 



V3 



f3' 



V6 



{D < 1). 



(15) 



We have also given the small D expansion which follows from 
(6). The plot of the exact function F = r(£') is shown in 
Fig. 2(b); notice that it is not D —D symmetric. 

In the lattice geometry of Fig. 1(b), the effective interaction 
between the pseudospins is determined by the spin-spin inter- 
actions, and appears already at order J'. From (12) we easily 
obtain 
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(ij> 



(16) 
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FIG. 3: (Color online) Types of magnetic order: (a) Chiral, (b) 
Coplanar, and (c) Coplanar in an external magnetic field (blue ar- 
rows are longer than green ones.) 



implying that ferromagnetic ordering takes place (since J' > 
0), i.e. on every tetrahedron (Tl^) = 1/2 (or (i;^) = 
— 1/2). The Ising symmetry T-^ —T^ is spontaneously 
broken. Physically, the state {T^) = 1/2 corresponds to 
a unique ground state formed as the linear combination: 
(1$) +i|5')) 1^/2. There is magnetic order in this ground 
state, as follows from (12): (Si) = (l/2)rAi; it is shown in 
Fig. 3(a). This is a magnetic "chiral" state, characterized by a 
non-zero scalar chirality (S^ • (Sj x Sfe)) 7^ 0. The magnitude 
of the magnetic moment is | (S^) | = (A/3/2)|r|, and the order 
is "weak" in a sense that F is determined by the value of the 
DM vectors. It is interesting to note that chiral order has been 
discussed in the context of spin liquid physics, where one 
can presumably have a state with broken time-reversal sym- 
metry and yet not magnetic order In our case however the 
presence of DM interactions (which break spin-rotation sym- 
metry) makes any time-reversal broken state magnetic. 

Magnetic chiral spin states were also discussed in the con- 
text of the full 3D pyrochlore lattice and in other (quasi 2D) 
non-frustrated situations. '^-^^ Such states generally compete 
with dimer order, which technically manifests itself in the 
presence of T^^T-^ operators in the effective Hamiltonian. 
Averages of such operators in the ground state lead to dimer 
order, and in turn diminish the magnetic, T^^ component. ^''^'^^ 
The specifics of this competition depend on the lattice. In 
our case (Fig. 1(b)) the issue of spontaneous dimerization has 
not been studied, to the best of our knowledge. However, it 
is clear that couplings involving T^, 2^ appear only in order 
(J')'', and higher Thus we will assume they can be neglected 
in the limit J' ^ 1. Given the coefficient in (16), a more 
precise criterion for the magnetic order to be dominant over 
potential spontaneous dimerization is J'F^ > {J')*, which 
we implicitly assume to be satisfied. Thus (16) determines 
the ground state structure within the degenerate subspace, and 
the degeneracy is lifted as explained previously, by locking 
{T^) — 1/2, which will be assumed from no on. (The choice 
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{Tf) = —1/2 leads to time-reversal of all magnetic states.) 

Now we write the effective Hamiltonian for the lowest ex- 
cited state h, Eq. (7). It is convenient to express the spin op- 
erators in the ground state {T^) = 1/2 via the 6, operators, 
similarly to the way it is summarized for the case of zero chi- 
rality {D = 0) in (Al). Performing the necessary calculations, 
we obtain 



•^1,3 



y2t&^ +h.c., = -Sl^. 



(17) 



In (17) we use notation such that the lower left index of 5*"^ 
(which stands for either Sf or Sf) corresponds to the up- 
per sign on the right hand side, and the lower right index — 
to the lower sign (i.e. Sf = (l - iVS) 6^ + h.c, and 

= (1 - iVs) W + h.c, etc.) The formulas (17) re- 
fer to a single tetrahedron which we label (as before) with a 
bold index i (i.e. h will carry this index, h bi.) The site 
index i = 1, 2, 3, 4 of the spins S'f on a tetrahedron again fol- 
lows the convention shown in Fig. 1(c). We also emphasize 
that equations (12),(17) simply mean that the spin operators 
have the same matrix elements as the right hand sides of those 
equations, i.e. the expressions should be added up to obtain 
the total spin operators. 

The coefficient t is defined as 



t ; 



P ^signP) 



6^2 2VG 6V2 



(18) 



and the small D limit expansion is also given. Figure 2(a) 
shows a plot of t. 

The effective Hamiltonian, describing the low-energy dy- 
namics of b can now be readily obtained: 

Hb = Aj^blh + J^i^Mbi+hblbl+h.c 



(ij> 



(19) 



The parameters appearing in Hi, are 

A = Eb-Eo, 

ti = -2t^{J' + V2D'), 

h = i(^)t{2V2J' + D'). 



(20) 



Here we have also included the DM interactions D' which ap- 
pear on inter-tetrahedral bonds (and we set D' = D from now 
on); their distribution is explained in Ref. 17. The presence 
of D' is not qualitatively (or even quantitatively) important 
for our following discussion. The quantity A is the on-site 
gap already discussed previously, while ti , t2 originate from 
the representation (17) for two spins on neighboring tetrahe- 
dra (thus ti, i2 ^ J't^-) Finally, ~ J'tT reflects terms of 
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FIG. 4: (Color online) (a) The h condensate density n = in 
zero field (solid red line), and in a finite magnetic field (dashed line), 
(b) The parameter F characterizing the level of coplanarity of the 
spin order at zero field (solid red line) and in a finite magnetic field 
along [001] (dashed line.) 



the type Ti{b\ +h.c.), i.e. the coupling between the magnetic 
component of the ground state pseudospin and the magnetic 
excited state (we have set TV {Tf) = 1/2, as (16) de- 
mands.) 

From (19) it follows that b condenses with a condensate 
value of 



- (6) 



^3 



A + 4(ti +t2)^ 



(21) 



calculated in the mean-field approximation. It is easy to check 
that (6) < 0. The b bosons are hard-core, in order to represent 
correctly the original spin operators. However their hard-core 
nature can be neglected, and thus the mean-field works well, 
as long as the condensate density n = is small. We 

find this to be the case for \D\ < 0.4; beyond this value the 
boson repulsion would stabilize the rise of (6), but this regime 
occurs only for unphysically large values of D. The density n 
is plotted in Fig. 4(a) with a solid red line, and from now on we 
set J' = J, for definitiveness. It is clear that a big difference 
exists between the two cases D > and D < 0, partially 
reflecting the difference in the behavior of A (Fig. 2(a)). For 
D > Q the value of (6) is significantly different from zero, 
while (6) « for L» < 0. 

Let us now investigate how the presence of (6) 7^ affects 
the spin order in the ground state. We combine equations (12) 
and (17), with = 1/2, b ^ (&), to obtain (the site index i 
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is defined as in Fig. 1(c)) 

(Sf> = {-iy+''^{l-F), z = l,2,3,4, 
{Sf) = = -{SI) = = + I 

(sf) = (SI) = = -{sy) = - £ (^1 + 1) 



,(22) 



where we have defined 



(23) 



The function F is plotted in Fig. 4(b)(solid red line.) For F = 
the magnetic order is of the chiral type (Fig. 3(a)), while for 
F = 1 it is coplanar, as shown in Fig. 3(b). At the coplanar 
point the magnetic moment is |(Si)| = (3-\/2/4)r. For an 
intermediate value of F the order in not coplanar We find the 
universal trend that for D < the order is almost chiral (F 
is small), and for D > there is a strong tendency towards 
coplanar order. 



B. Enhancement of coplanar order by magnetic field 

We now proceed to investigate how an external magnetic 
field can affect the ordering tendencies described above. It is 
known, for example in dimer systems, that a field can in- 
duce magnetization perpendicular to it in the presence of DM 
interactions. The exact solution of the problem in our case is 
rather complex, and below we only give a summary of the re- 
sults for the case of weak fields. We consider a magnetic field 
H along [001], i.e. in the z direction in the coordinate system 
of Fig. 1(c). The field is assumed to be weak in the sense that 
H ^ D. The presence of a field leads to changes in the spec- 
trum, and we have found that quantitatively the most visible 
effect is related to the change of the h level energy. In fact in 
magnetic field the b state mixes with the state \Pz) (Eq. (9)), 
and the new eigenstate, ci\b) + C2\Pz), has energy 
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y 



-,1/2 



(24) 



The on-site gap in a field, /S.h ~ Ei,{H) — Eq, is plotted in 
Fig. 2(a) with a dashed line (for the specific choice H ~ 2D). 
This behavior in turn enhances the tendency towards coplanar 
order, as evidenced by the dashed lines in Fig. 4(a,b). We do 
not present the exact form of the spin operators in this case, 
but just mention that they differ slightly from the zero-field 
case (22). In particular, at the coplanar point, the magnetic 
moments are are not equal on all sites of the tetrahedron, as 
shown in Fig. 3(c), where the blue (upper) arrows are longer 
than the green (lower) ones. 



two main types of spin order are possible, as summarized in 
Figure 3. These are (non-coplanar) chiral and coplanar order- 
ings, depending on the sign of D. In our model the spin or- 
der also generally deviates slightly from the exact chiral and 
coplanar configurations, and this deviation itself depends on 
D. It is interesting to note that magnetic order of the kind de- 
scribed in this work was also found in the Sp{N) (large N) 
approach to the Heisenberg pyrochlore model. In addition, 
collinear order is also possible in that case, while DM interac- 
tions do not favor collinear states. The types of DM-induced 
order (chiral and coplanar) we found in the quantum case are 
also consistent with the Monte-Carlo results for the classical 
(3D) model. 1^ 

It is also important to emphasize that our approach assumes 
a specific type of lattice-symmetry breaking (as the explicit 
anisotropy of Fig. 1(b) demands.) In the full (3D) quantum 
pyrochlore lattice, with Heisenberg interactions only, it is be- 
lieved that (spontaneous) lattice-symmetry breaking always 
takes place in the ground state, and a certain dimeriza- 
tion pattern sets in. Two-dimensional projections previously 
studied"''^ also exhibit valence-bond solid ground states. In 
our quasi-2D version (Fig. 1(b)), the dimerization tendency 
on the tetrahedra is very weak (as it occurs only in fourth 
order of perturbation theory) and does not interfere with the 
DM-induced magnetic order formation. At the same time our 
results are not directly relevant to the 3D pyrochlore lattice. 
Also, in this work we have not provided a dynamical mech- 
anism for the layer decoupling, which we have taken as our 
starting point (the presence of spin-phonon interactions with 
the correct symmetry can certainly accomplish this task.) In- 
stead, our goal has been the study of time-reversal symmetry 
breaking in certain anisotropic model situations with strong 
frustration. 

Finally, the magnetic order sets in only below the Ising tran- 
sition temperature Tc ^ JT^, as dictated by the pseudospin 
interactions. Since F is determined by the DM interaction 
(Fig. 2(b)), we expect Tc to be small, and so are the magnetic 
moments (in all magnetic patterns), |(Si)| ~ |r|. So far we 
are not aware of any convincing evidence that DM-induced 
order takes place in the pyrochlore -related compounds; nev- 
ertheless the presence of genuine antiferromagnetic order is 
a fundamental property of pyrochlore systems, which distin- 
guishes them from other situations with more "trivial" mani- 
festations of DM interactions (such as weak ferromagnetism.) 
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IV. CONCLUSIONS 

In this work we showed that, in the anisotropic version of 
the pyrochlore lattice and in the presence of DM interactions. 



APPENDIX A: SPIN OPERATORS ON A TETRAHEDRON 

For completeness we summarize the values of the various 
matrix elements used in the main text. AtD = 0, |si), \s2) are 
the singlet ground states, and \ta), \pa), \qa), = x,y, z are 
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S = 1 states. The lower left (site) index on the spin operators 
corresponds to the upper sign on the right hand side (and the 
lower right index corresponds to the lower sign, if different.) 

= -^tisi ± ^^pUi T ^gis2 +h.c. 

(Al) 

The states are defined as 



\s^) = -= [| Tiit) + I itti) + I TTii) + I iitt) 
-2(1 TiTi) + UTiT»], 



\S2) 




TiiT) + 


liTTi)- 


1 mi) - 


1 ilTT)], 


\t.) 
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[| nn>- 


-unt)] 










Tiit) - 


1 iTTi> + 


1 TTii) - 


1 iitt)], 






-1 TiiT)- 


f litn)- 


f 1 tni)- 


-UiTT)] 



(A2) 



Only the spin zero components of the triplets are shown, and 
in the notation of the type | tiit) the spins follow the order 
1, 2, 3, 4, with the site labels defined in Fig. 1(c). 
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